D — D mixing in the framework of the HQE revisited 
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We reconsider the leading HQE contributions to the absorptive part of the mixing amphtude of neutral D 
mesons by taking also as corrections and subleading 1/mc corrections into account. We show that these 
contributions to do not vanish in the exact SU(3)f limit and they also can have a large phase. Moreover, 
we give an example of a new physics model that can enhance these leading HQE terms to ri2, which are orders 
of magnitude lower than the current experimental expectation, by a factor in the upper double-digit range. 
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Introduction. — Mixing of neutral mesons provides an 
excellent testing ground for the standard model (SM) and its 
possible extensions. In the system the following quantities 
have been measured [ 1 ] 

Ar 7.3 ±1.8 A7\/ 9.1^2^ 
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Mixing formalism. — The mixing of neutral mesons is 
described by box diagrams with the absorptive part ri2 and 
the dispersive part A/12. The observable mass and decay rate 
differences are given by (0 :— arg[— A/i2/ri2]) 

{AMf~UATf = 4|A/i2|2 



4 

AA/Ar 



|ri2p, 



(2) 



4|A/i2||ri2|cos(0). 

If |ri2/Afi2| < 1, as in the case of the i?s system (« 5-10"^) 
or if ^ 1, one gets the famous approximate formulae 

AA// = 2|A/i2|, Ar = 2|ri2|cos(/). (3) 

The experimental values for x and y suggest that in the Z?" 
system |ri2/A/i2| « the size of the mixing phase cf> 

will be discussed below. 

Leading HQE predictions. — The absorptive part of the 
box diagram with internal s and d quarks can be decomposed 
according to the CKM structure as 

ri2 = - (A^r,, + 2A,Adr,d + xITm) , (4) 

with Xx = VcxV*^. The application of the heavy quark ex- 
pansion (HQE), which turned out to be very successful in the 
B system, to the charm system typically meets major doubts. 
Our strategy in this work is the following; instead of trying 
to clarify the convergence of the HQE in the charm system in 
advance, we simply start with the leading term and determine 
corrections to it. The size of these corrections will give us 
an estimate for the convergence of the HQE in the _D" system. 
To this end we first investigate the contribution of dimension-6 
(D=6) operators to ri2, see Fig. (1). Next we include NLO- 
QCD corrections, which were calculated for the Bs system 
[2, 3, 4, 5] and subleading terms in the HQE (dimension-7 
operators), which were obtained in [6, 7]. To investigate the 
size of as and l/rric corrections in more detail we decom- 
pose Tss.sd.dd into the Wilson coefficients G and Gs of the 
AC = 2 operators Q and Qs (for more details see [5]) 
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FIG. 1: Contributions to from operators of dimension 6 (D=6). 
The leading order QCD diagram is shown in the left panel, an exam- 
ple for Os corrections is shown in the right panel. 



with Q = M„7'^(l - 75)Ca ■ U;37m(1 - l5)C(3, 

Qs = Wa(l + 75)Ca ■ + 75)c/3. 



(6) 



The effect of the QCD corrections has already been discussed 
in [8]. In our numerics we carefully expand in ag'. the leading 
order QCD contribution consists of leading order AC = 1 
Wilson coefficients inserted in the left diagram of Fig. (1), 
while our NLO result consists of NLO AC = 1 Wilson coef- 
ficients inserted in both diagrams of Fig. (1) and consistently 
throwing away all terms which are explicitly of 0{a^). Fol- 
lowing [9] we have also summed terms like 2 In z to all orders, 
therefore we use in our numerics z = fhs{fhc)'^ /mcifricY ~ 
0.0092. The matrix elements in Eq. (6) are parameterized as 



{D\Q\D) 
{D\Qs\D) 



-fnMlBif,) 
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We take fo = 212(14) MeV from [10] and we derive 
B{nic) = 0.9±0.1 and Bs(mc) ~ 0.9 from [11, 12]. Finally 
we use the A/5 scheme for the charm mass, fhc{rhc) = 1.27 
GeV. For clarity we show in the following only the results for 
the central values of our parameters, the error estimates will 
be presented at the end of the next section. We obtain 
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For the error estimate we vary fii between 1 GeV and 2mc 
Combining G and Gs to T^x we get (in units of ps~^) 





LO 


ANLO-QCD 


Al/mc 


E 




2.68 


-0.68(-25%) 


-0.76(-28%) 


1.24 




2.70 


-0.67(-25%) 


-0.76(-28%) 


1.26 




2.71 


-0.66(-24%) 


-0.76(-28%) 


1.29 



Using instead the operator basis suggested in [5] with 

Qs = ■Sq(1 +75)c/3 • 11/3(1 + 75)Cq , 



(8) 



{D\Qs\D} = blMh ( , .^i"" , , 
and Bs{mc) ~ 1.1 leads to 





LO 


ANLO-QCD 


Al/mc 


E 




1.71 


+0.02(+l%) 


-0.34(-20%) 


1.40 




1.73 


+0.04(+2%) 


-0.34(-20%) 


1.42 




1.74 


+0.05(+3%) 


-0.34(-20%) 


1.45 



All in all we get large QCD (up to 50%) and large llrric cor- 
rections (up to 30%) to the leading D=6 term, which consider- 
ably lower the LO values. Despite these large corrections, the 
HQE seems not to be completely off. From our above inves- 
tigations we see no hints for a breakdown of OPE. The same 
argument can be obtained from the comparison of B and D 
meson lifetimes. In the HQE one obtains 
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where the leading term Fq cx c^ckm corresponds to the 
free quark decay and all higher terms in the HQE are com- 
prised in sr. For the ratio Fo,D/ro.B one gets a value close 
to one. Higher order HQE corrections in the B system are 
known to be smaller than 10 % [13]. Using the experimental 
values for the lifetimes we get 



— « 1.4.. .4 (Exp.) w 1 
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(10) 



From this rough estimate one expects higher order HQE cor- 
rections in the D system of up to 300 %. So clearly no pre- 
cision determination will be possible within the HQE, but the 
estimates should still be within the right order of magnitude. 

Cancellations. — As is well known huge GIM cancella- 
tions [14] arise in the leading HQE terms for D mixing. To 
make these effects more obvious, we use the unitarity of the 
CKM matrix (A^ + As + A;, = 0) to rewrite the expression for 
the absorptive part in Eq. (4) as 

(11) 

Note that the CKM structures differ enormously in their nu- 
merical values: Xd.s oc A and At oc A^ [15] in terms of the 
Wolfenstein parameter A « 0.2255. In the limit of exact 
SU{3)f symmetry, Tgs — ^sd — ^dd holds and therefore. 



contrary to many statements in the literature, F12 — — 
is not zero but strongly CKM suppressed. Next we expand the 
arising terms in z. Using Eq. (5) we get in LO 



^ So 

rD=6 
^ sd 



2.71252- 3. 14543Z 
2.71252- 1.57271Z 



9.98427r 
4.99214^2 



(12) 



The first term in the above equations obviously corresponds 
to F^j^^. For the combinations in Eq. (1 1) we get 

^D=6 ^r. ir,r^^-3 ^ \7.3 



(r. 



2rsd 

[Tsd 



+ Tdd) 

~ ^dd) 



D=6 



= -26.1794r 
= -1.57271z; 



A 



A' 
2.8 



(13) 



To make the comparison with the arising CKM structures 
more obvious, we have expressed the size of these combina- 
tions also in terms of powers of the Wolfenstein parameter A. 
As it is well known, we find in the first term of Eq. (11) an 
extremly effective GIM cancellation, only terms of order 
survive. In NLO we get 



= 1.2903 - 5. 35346Z- 



.77624z' 



D=6J 
sd 



= 1.2903 - 2. 67673z-4.68377r 



The arising combinations in Eq. (11) read now 



D=6,7 



0.59r 



{Tsd-rddf-'''' = -2.68Z 



43.6r 
« A^ 



+ .... (14) 

X' ' - A^-s, 
(15) 



The fact that now the first term of Eq. ( 1 1) is of order com- 
pared to z^ in the case of the LO-QCD value was discussed in 
detail in [ ] and later on confirmed in [8]. These numbers are 
now combined with CKM structures, whose exact values read 



Xd 

Xs 

with c 



-C12C23C13S12 



C?2Cl3S23Sl3e"^' 



+ C12C23C13S12 - Si2Cl3S23Sl3e 

ci3S23Si3e''^^ =0(A5 + zA^) 



o{x' 



(16) 



cos 



and Sij ~ sin( 



). Looking at Eq. (16), it 
is of course tempting to throw away the small imaginary parts 
of Xd and A,, but we will show below that this is not justified. 
Doing so and keeping only the leading term in the CKM struc- 
ture (C12C23C13S12), which is equivalent to approximate Af, = 
0, one gets a real ri2 which vanishes in the exact SU (3) p 
limit. Keeping the exact expressions, we see that the first term 
in Eq. (1 1) is leading in CKM (O [A^ + iX^] ) and has a negU- 
gible imaginary part, but it is suppressed by 0.6z^ ~ A^ ''. The 
second term in Eq. (11) is subleading in CKM (O [X^ + iX^] 
and it can have a sizeable phase. This term is less suppressed 
by SU{S)f breaking (« 2.7z w A^-^). The third term in 
Eq. (11) is not suppressed at all by SU{3)f breaking, but it 
is strongly CKM suppressed {O [A^" + iX^^] ). For clarity we 
compare the different contributions of Eq. (1 1) in the follow- 
ing table 





1st term 


2nd term 


3rd term 


LO 


26A^z^ 
« A9-3 


3.1AsAbZ 
« A8-4 


2.7A^ 
« A9-3 


NLO 


0.59A^z^ 
« A8-7 


5.4:XsXbZ 

« A^ o 


0.66A^ 

« Ai"-3 
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From this simple power counting, we see that a priori no con- 
tribution to Eq. (11) can be neglected. Having the comment 
[15] in mind, we find however, that the first two terms of 
Eq. (11) are of similar size, while the third term is suppressed. 
Moreover, the second term can give rise to a large phase in 
ri2, while the first term has only a negligible phase. To make 
our arguments more solid we perform the full numerics using 
the CKM values from [ I 7] and obtain for the three contribu- 
tions of Eq. (11) 

10^rf2=^'^ = -0.38 + 0.00002i (1st term) 

-6.34 - 15.0i (2nd term) 

+0.19 - O.m (3rd term) 

= -6.53 - 15.2i (11...38)e-*(o-5 -2-3). 

Here we show for the first time the errors, they are estimated 
by varying /^ii between 1 GeV and 2toc and the dependence 
on the choice of the operator basis. The first term in Eq. (11) 
turns out to be very sensitive with respect to the the exact val- 
ues of the bag parameters and its real part can, within errors, 
be numerically of the same size as the second term, which 
features a large imaginary part. Furthermore, even the third 
term can give a non-negligible contribution, in particular to 
the imaginary part. 

To summarize, we have demonstrated that the typical ap- 
proximation Ab « 0, which is equivalent to neglecting the 
imaginary parts of \d and is wrong for the case of the lead- 
ing HQE prediction for and yields the wrong conclusion 
that cannot have a sizeable phase. We get for the first 
terms in the OPE a value for hd of 

These values are still a factor of 0.8. ..1.3 ■ 10'' smaller than 
the experimental number. This is in contrast to our previous 
expectations that the HQE should give at least the right order 
of magnitude. Moreover, we do not confirm the observation 
made in [S] that the NLO result for is almost an order of 
magnitude larger than the LO result. 

Higher HQE predictions. — In [18, 19, 20] higher or- 
der terms in the HQE of D mixing were discussed. If the 
GIM cancellation is not as effective as in the leading HQE 
term, operators of dimension 9 and dimension 12, see Fig. (2), 
might be numerically dominant. In order to obtain an imag- 
inary part of the loop integral, the operators of dimension 9 
have to be dressed with at least one gluon and the operators 
of dimension 12 with at least two gluons. If we normalize the 
leading term (left figure of Fig. (1)) to 1, we expect the D=9 
diagram of Fig. (2) to be of the order 0{as{4:7T){qq) /m^) « 
0.03 and the D=12 diagram of Fig. (2) to be of the order 
C'(a^(47r)^(^g)^/m^) « 10^^. As explained above, the for- 
mally leading term of D=6 is strongly GIM suppressed to a 
value of about 2-10-^ and the big question is now how severe 
are the GIM cancellations in the D=9,12 contribution. For the 



{dd),{ss) {dd),{ss) 




d,s {dd),{ss) 



FIG. 2: Contributions to from operators of dimension 9 (D=9, 
left panel) and dimension 12 (D=12, right panel). To obtain an imag- 
inary part the D=9 diagrams have to be dressed with at least one 
gluon and the D=12 diagrams with at least 2 gluons. 



contributions to y^i we get the naive expectations 
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If there would be no GIM cancellations in the higher OPE 
terms, then the D=9 or D=12 contributions could be orders 
of magnitudes larger than the D=6 term, but in order to ex- 
plain the experimental number still an additional numerical 
enhancement factor of about 15 has to be present. For more 
substantiated statements r{^~^'^^ has to be determined explic- 
itly, which is beyond the scope of this work [2 1 ]. This calcu- 
lation is also necessary in order to clarify to what extent the 
large phase in ri2 from the first OPE will survive. In order to 
determine the possible SM ranges of the physical phase (f), in 
addition one has to determine M12. 

New physics. — Finally we would like to address the 
question, whether new physics (NP) can enhance ri2. In the 
Bs system it is argued [22] that ri2 is due to real intermediate 
states, so one cannot have sizeable NP contributions. More- 
over, the mixing phase in the Bg is close to zero, so the cosine 
in Eq. (3) is close to one and therefore NP can at most modify 
(f), which results in lowering the value of AF compared to the 
SM prediction. In principle there is a loophole in the above 
argument. To ri2 also AS = 1 penguin operators contribute, 
whose Wilson coefficients might be modified by NP effects. 
But these effects would also change all tree-level B decays. 
Since this is not observed at a significant scale, it is safe to say 
that within the hadronic uncertainties ri2 = ^12^ there- 
fore the argument of [ J :] holds. Since in the system the 
QCD uncertainties are much larger, also the possible effects 
might be larger but not dramatic. The peculiarity of the 
system - the leading term in the HQE is strongly suppressed 
due to GIM cancellation - gives us however a possibility to 
enhance ri2 by a large factor, if we manage to soften the GIM 
cancellation. This might be accomodated either by weakening 
the SU(3) F suppression in the first two terms of Eq. (11), see, 
e.g., Petrov et al. [23], or by enhancing the CKM factors of the 
last two terms inEq. (11). The latter can be realized in a model 
with an additional fourth fermion family (SM4). The usual 
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FIG. 3: The enhancement factor /TiT using the possible val- 
ues for VcKMi found in [16]. The color encoded scale denotes the 
logarithm of the number of allowed parameter points of Vckam- 



CKM matrix is replaced by a four dimensional one (Vck am) 
and the unitary condition now reads Xd + Xs + + Xf = 0. 
Eq. (11) is replaced by 

ri2 = ^xlir,,^2T,d + Tdd) (18) 

+2As(Ab + Xb') (Tsd — Tdd) — (Ab + Xi,/)^Tdd ■ 

In [ 1 6] an exploratory study of the allowed parameter space of 
Vck MA was performed and as expected only very small mod- 
ifications of Xd and are experimentally allowed. In almost 
all physical applications these modifications are numerically 
much smaller than the corresponding hadronic uncertainties 
and therefore invisible. However, in the D'-* mixing system it 
might happen, that all dominant contributions cancel and only 
these modifications survive. In the SM the first two terms of 
Eq.(ll) are numerical equal. In the SM4 the numerical hi- 
erarchy depends on the possible size of Xb', see Eq.(18). In 
particular, it was found in [16] that currently a value of Xb' 
of the order A"^ is not excluded. This means that the second 
term of Eq.(18) could be greatly enhanced by the existence of 
a fourth family and also the third term would now become rel- 
evant. Using experimentally allowed data points for Vckam 
from [16] we have determined the possible values of in 
the SM4 from Eq. (18); enhancement factors of a few tens are 
easily possible, see Fig. (3). 

Discussion and Conclusions. — In this work we inves- 
tigated the leading HQE contribution to the absorptive part 
of D mixing and corrections to it. We found that the size 
of these corrections is large, but not dramatic (« 50% QCD, 
« 30% 1/toc). So it seems that the HQE might be appropriate 
to estimate the order of magnitude of ri2. We have explained 



that r^^"^'*^ gives, due to huge GIM cancellations, a value of 
y which is about a factor of about 10000 smaller than the cur- 
rent experimental expectation, but it can have a large phase 
and it also does not vanish in the exact SU(3)i? limit. Due 
to this peculiarity, it might be possible that the HQE result is 
dominated by D=9 and D=12 contributions, if there the GIM 
cancellations are less pronounced, but to make more profound 
statements - in particular about the standard model value of 
the mixing phase - these higher dimensional coiTections have 
to be determined explicitly. Finally we have shown that new 
physics can enhance T12 by a high double-digit factor. 
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